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IN  SITU  CHARACTERIZATION  OF  SOILS  FOR  PREDICTION 
OF  STRESS-STRAIN  RELATIONSHIP 

1.  Introduction 

The  solution  of  most  geotechnical  engineering  problems  requires  a 
knowledge  of  engineering  properties.  One  of  the  major  difficulties  in  obtaining 
soil  properties  accurately  is  the  disturbance  during  sampling  or  the  use  of 
penetration  methods  of  in  situ  testing  which  alters  the  engineering  properties. 
In  si*u  testing  techniques,  however,  have  begun  to  play  an  increasingly  important 
role  in  tine  determination  of  soil  properties.  This  trend  towards  in  situ  testing 
techniques  could  be  ascribed  to  several  reasons  as  follows:  (1)  the  soil  is  tested 
in  its  in  situ  environmental  conditions  which  influence  the  Engineering  properties, 
(2)  continuous  data  through  the  profile  may  be  obtained,,  and  (3)  properties  can 
be  obtained  in  cases  where  obtaining  undisturbed  samples  is  very  difficult,  such 
as  in  the  case  of  saturated  sands. 

Presently,  several  in  situ  testing  techniques  are  employed  for  the  evaluation 
of  soil  properties  and  for  the  establishment  of  empirical  criteria  for  the  prediction 
of  potential  behavior  of  soils.  The  standard  penetration  test  (SPT)  and  the  cone 
penetration  test  (CPT)  are  the  most  widely  used  in  situ  testing  techniques. 
Other  in  situ  testing  techniques  such  as  the  pressure  meter  test  (PMT),  the  Iowa 
Bore  Hole  Shear  Test  (BST)  and  the  dilatometer  are  being  used  on  a  modest 
scale. 

The  reliability  and  usefulness  of  test  results  obtained  from  the  in  situ 
testing  techniques  described  above  are  limited  owing  to  various  reasons.  The 
drainage  conditions  in  soils  during  the  in  situ  testing  may  be  unknown  in  some 
tests.  The  failure  modes  may  not  simulate  those  anticipated  for  the  actual 
project  and  in  some  cases  the  exact  failure  mode  is  unknown.  Due  to  these 
unknown  drainage  and  failure  conditions,  and  the  specific  nature  of  testing 
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procedures  and  the  failure  mode  simulated  in  different  techniques,  interpretation 
of  the  test  results  is  highly  empirical  at  the  present  time.  The  generalization 
of  the  test  results,  obtained  from  the  current  in  situ  testing  techniques  in  order 
to  analyze  soil  behavior  under  general  loading  conditions  and  various  drainage 
and  boundary  conditions,  is  very  difficult: 

A  nondestructive  method  of  chaiacterizing  particulate  systems  is  presented 
by  considering  the  electrical  properties  of  soils  which  can  be  determined  in  situ 
without  causing  disturbance  to  the  soils.  The  objectives  of  this  report  are: 

1)  characterization  of  soil  by  electrical  methods 

2)  quantification  of  the  compositional  and  aggregate  properties  of  soils 
in  terms  of  electrical  parameters 

3)  quantification  of  the  inter-  and  intra-cluster  Void  ratios 

4)  establishment  of  correlations  between  electrical  parameters  and 
compression  index,  X,  swell  index  K,  and  the  slope  of  the  critical 
state  line,  M  for  soils 

5)  methodology  to  predict  stress  strain  relationship  of  soils  in  situ 

6)  establish  correlation  between  KQ  and  the  electrical  parameters. 

2.  Characterization  of  Soils  bv  Electrical  Methods 

When  an  alternating  electrical  field  is  applied  to  a  clay-water -electrolyte 
system,  a  response  is  produced  which  can  be  measured  in  terms  of  a  resistance, 
R,  and  a  capacitance,  C.  The  measured  value  of  the  capacitance  can  be 
converted  into  a  quantity  known  as  the  dielectric  constant.  This  value  is  defined 
as  C/^.Q,  in  which  Co  is  the  capacitance  of  a  condenser  with  oniv  a  vacuum 
between  the  electrodes.  The  dielectric  constant  is  actually  a  measure  of  the 
ability  of  the  clay  to  store  electrical  potential  energy  under  the  influence  of 
an  electric  field.  From  a  knowledge  of  the  dimensions  of  the  sample  the 


dielectric  constant  and  conductivity  can  be  calculated  from  values  of  R  and  C 
using  the  following  relationships. 


(1) 


(2) 


in  which  d  =  length  of  a  specimen;  A  =  cross-sectional  area;  and  e  =  the 

a 
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dielectric  constant  oi  vacuum  (8.85  x  10"  farad/cm). 

The  dielectric  constant  of  a  dry  silicate  material  is  about  4.5,  and  that 
of  water  is  about  79.  A  mixture  of  soil  and  water  should,  therefore,  have  a 
dielectric  constant  weighted  between  4.5  and  79.  However,  when  the  dielectric 
constant  of  a  clay-water-electrolyte  system  is  measured  with  an  alternating 
current  in  the  radio-frequency  range,  it  is  found  to  be  in  excess  of  the  sum  of 
the  weighted  dielectric  constant  of  the  components.  This  measured  value,  e', 
referred  to  as  the  "apparent  dielectric  constant,"  reflects  the  heterogeneous 
nature  of  the  path  of  the  current  and  the  electrical’ properties  of  the  pore  fluid 
and  the  clay  mineral  as  shown  by  Arulanandan  and  Smith  (1973). 

When  the  conductivity,  c,  and  apparent  dielectric  constant,  e',  of  a 
cohesive  system  are  measured  as  a  function  of  frequency,  in  the  radio  frequency 
range,  it  is  found'  that  e'  and  a  vary  as  shown  in  Figure  2.  This  variation  of 
o  and  e'  is  referred  to  as  electrical  dispersion.  However,  for  a  granular  system 
it  has  been  observed  that  a  and  £'  are  independent  of  frequency  (see  Figure  1). 
It  is  thus  possible  to  classify  all  particulate  systems  into  two  groups.  All 
cohesive  systems  are  characterized  by  electrical  dispersion  behavior  and  granular 
systems  by  nondispersive  electrical  behavior. 


This  method  of  classification  is  non-destructive,  quantitative  and 
fundamental  in  nature  as  opposed  to  the  present  qualitative,  arbitrary  methods 


of  classifying  soils  based  on  grain  size  and  plasticity  indices. 

The  factors  that  influence  the  magnitude  of  dielectric  dispersion  A eQ  (see 
Figure  2)  were  investigated  in  detail  by  Arulanandan  et  al.  (1973).  It  was  shown 
that  AeQ  depend  predominately  on  the  type  and  amount  of  clay  mineral  and  the 
fabric  having  a  second  order  effect.  Thus  the  composition  of  cohesive  soils  can 
be  characterized  nondestructive^  and  quantified  in  terms  of  the  parameters 
A  eQ.  The  significance  oi  this  parameter  to  the  classification  of  soils,  Arulanandan 
et  al.  (1973),  to  the  prediction  of  cation  exchange  capacity  Fernando  et  al. 
(1975),  swell  potential  Basu  and  Arulanandan  (1973),  lime  reactivity  Castle  and 
Arulanandan  (1979),  erosion  potential  Heinzen  and  Arulanandan  (1977)  and 
compression  index  Scharlin  (1974),  has  been  demonstrated. 

2.1  Electrical  Models  to  Predict  Electrical  Dispersive  Behavior 

Arulanandan  and  Mitra  (1970),  Arulanandan  and  Smith  (1973)  have  shown 
it  is  possible  to  describe  the  electrical  dispersion  behavior  of  soils  in  terms  of 
an  electrical  network  model. 

The  current  through  the  sample  is  considered  to  have  three  paths  (Figure  3): 
(1)  through  the  intercluster  solution  and  clusters,  (2)  through  clusters  in  contact 
with  each  other,  and  (3)  through  the  intercluster  solution.  The  corresponding 
three  element  model  representation  of  the  impedances  of  individual  components 
is  shown  in  Figure  4c. 

Referring  to  Figure  4c,  kp  =  conductivity  of  the  solution,  ef  =  dielectric 
constant  of  the  cluster,  e$  =  dielectric  constant  of  the  solution,  a  =  area  of 
cluster -solution  path,  b  =  area  of  the  particle  contact,  path  and  c  =  area  of  the 
solution  path.  In  the  case  of  soils  saturated  in,  water,  e$  is  approximately  79 
and  ks  is  approximately  equal  to  the  conductivity  of  the  pore  fluid. 


The  expressions  for  the  apparent  dielectric  constant,  e',  and  the  apparent 
conductivity,  o’,  have  been  .jui  u  utilizing  the  electrical  network  theory  as 
follows: 

2  2  2  2 
ek’  ek  __ee.ee 
a  r  r  s  s  r  2J2  ,  r  s  r  s^  __ 

£  =  dTTdTs  {rr  +  —  +  u  a  (u r  +  — +  ber +  ces  (3) 


kk2  k2k  _  e2k  eV 
a  r  r  s  r  s  J.  ,  s  r  r  s*i  Ll. 

dTudTs  +  —  +  wa  +  —V  +  bkr +  d<s 


where 


l>  3  2IIf 


-  ,ks  kr x2  22  £s  er\: 

s  =  (ud  +  d~]  +  “  a  +  T] 


f  =  current  frequency 

a  =  capacitance  of  a  unit  capacitor  in  vacuum  (8.85  x  10"^  farad/cm) 

kr  =  conductivity  of  the  cluster 

er  =  dielectric  constant  of  the  cluster 

k  s  conductivity  of  the  solution 
s 

e  =  dielectric  constant  of  the  solution 
s 

a  =  x-sectional  area  of  cluster-solution  path 

b  =  v-sectional  area  of  particle  contact  path 

c.  =  x-sectional  area  of  solution  path 

d  =  length  of  path  in  the  solids  in  the  solid-solution  phase. 

2.2  Quantification  of  Inter-  and  Intra-cluster  Void  Ratios 

Anandarajah  (1982),  Arulanandan  et  al.  (1983)  equated  the  three  element 
electrical  model,  Figure  4a,  described  in  the  previous  pages  to  the  model  shown 
in  Figure  4b  in  order  to  obtain  the  volume  properties  such  as  porosity. 


•  •  •  •  -  •  m  •  V  •  '  «  ,  ■  ,  •  .  '  .  • 


Defining; 


cluster  tortuosity,  t  j 

L2e  Llle 
"  L  'Lll 

(5) 

solution  tortuosity,  t^ 

L3e  L12e 
’  L  'L12 

(6) 

A1 

(7) 

3  "  AV2 

h  ^ 

b  =  ““ 2 
At? 

■ 

(8) 

*3 

c  =  — 

At2 

■ 

(9) 

but  A.  s  Aj  +  Aj 

+  Aj,  i.e 

.  atjt2  +  bt*  +  ct2  =1 

(10) 

A  “'He 

d  = 

,  1  - 

.  L12e 

d  =  i 

2 

(11) 

Therefore 

. 

L11  =  L’d 

and 

=  (I-d)L 

(12) 

and  assuming  electrical  volumes  equal  to  true  volumes. 

Then 

■V  •  L2e  2 

l  •  A2 

(13) 

II 

X 

/ 

L  >  Aj 

V*3  *  U<2 

(i») 

Alle  *  Llle  =  L11 

*  A1 

(15) 

A12e  *  L12e  s  L12 

*  A1 

A12e/At  1  1/,2 

(16) 
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Therefore,  assuming  e$  =  79  and  =  4.5 

gr  -  4»5 
nI  “  74.5 

but 


Hence, 

ef  -  4.5 
eI  =  79~ 


(21) 


(22) 


(23) 


and 


a  _VE_  =  CP 

Vp*Vc  l+ep>e j 


(l+ej) 


(24) 


(25) 


but 


LA 


(26) 


np  =  a( l-d)t jtj^ctj  (27) 

A  digital  computer  optimization  program  utilizing  the  simplex  method  was 
employed  to  analyze  the  dielectric  constant  and  conductivity  dispersion  curves. 
The  optimization  program  determines  the  values  of  the  geometrical  (a,  b,  c,  d), 
compositional  (€f,  kr,  k^,  zj  and  tortuosity  (tj,  t J  parameters  of  the  model 
shown  in  Figure  4c  such  that  the  frequency  dispersion  curves  of  the  apparent 
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dielectric  corstant  and  conductivity  calculated  from  the  theoretical  Equations  17 
and  13  fit  the  the  experimental  results  within  acceptable  error  limits,  the 
dielectric  constant  of  water,  was  kept  constant  at  79  during  ail  computer  rims. 

There  are,  in  fact,  eignt  random  parameters  b,  c,  d,  er,  kr,  k^,  tj,  t2 
involved  in  the  curve  fitting  procedure  by  optimization  technique.  At  first  sight, 
it  might  seem  that  the  choice  of  eight  seemingly  random  parameters  would  be 
enough  to  fit  almost  any  set  of  curves  and  that  the  physical  meaning  of  the 
result  is  doubtful.  It  is  important  to  realize,  however,  that  only  the  three 
independent  geometrical  parameters  are  genuine  unknowns  (b,  c,  d)  and  as  b  *  0,' 
there  are  only  two  independent  unknown  parameters,  "c"  and  "d".  The  respective 
conductivities,  kr  and  k$,  of  the  cluster  and  the  interstitial  solution  and  the 
dielectric  constant  of  the  solid  particles  must  have  values  that  are  reasonable 
from  the  point  of  view  of  physical  science.  One  would  expect  the  conductivity 
of  the  interstitial  solution  to  be  of  the  same  o.der  as.  the  conductivity  of  the 
pore  fluid  extract;  the  conductivity  of  the  solid  should  be  of  the  same  order  as 
the  isoconductivity  value  (Arulanartdan  et  al„  1973);  the  dielectric  constant  of 
the  clusters  should  be  between  4.5  (the  dielectric  constant  of  dry  silicate  mineral) 
and  79  (the  dielectric  constant  of  pure  water),  depending  on  the  water  content.' 
Although  these  limitations  were  not  written  into  the  computer  program,  the 
numerous  results  obtained  in  earlier  investigations  show  that  the  requirements 
were  indeed  met.  However,  a  large  number  of  computer  runs  is  necessary  to 
select  the  best  set  of  parameters  that  fit  the  experimental  curves  and  are 
reasonable  from  a  physical  science  point  of  view.  Thus  the  three-element' 
electrical  mode)  can  be  used  to  quantify  the  inter-  and  intra-cluster  void  ratios. 

The  quantification  of  the  inter-  and  intra-cluster  void  ratios  by  the  above 
approach  is  time-consuming,  as  several  computer  runs  are  needed  h  be  made 
to  obtain  the  model  parameters,  combined  with  a  certain  amount  of  judgement 
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needed  to  choose  the  right  parameters.  This  is  due  to  the  fact  that  Equations  17 
and  IS  are  highly  non-linear.  In  order  to  make  the  optimization  technique  more 
efficient  and  to  choose  the  suitable  model  parameters,  it  was  thought  necessary 
to  impose  an  additional  constraint  in  the  optimization.  This  would  be  achieved 
utilizing  the  measured  water  content  of  the  soil.  Once  the  dielectric  Constant 
and  conductivity  of  the  sample  in  the  radio  frequency  range  is  available,  to  get 
the  model  parameters  and  hence  the  inter-  and  intra-cluster  void  ratios  one 
needs  to  optimize  the  results  due  to  the  nonlinearity  of  the  problem. 

For  this  purpose  optimization  program  called  "Simplex"  is  used  with 
modification  to  suit  the  problem. 

Subroutine  Simplex  is  called  from  main  program  "Opti"  and  cos'*  function 
"Funct"  is  called  by  Simplex. 

In  the  cost  function  "Funct"  using  the  predicted  model  parameters 
dielectric  constant  and  conductivity  are  calculated  in  the  frequency  range  and 
compared  with  measured  to  obtain  the  minimum  percentage  error.  To  this  cost 
function  another  term  is  added  to  take  into  account  the  difference  between 
calculated  and  measured  water  contents. 

There  are  several  constraints  imposed  on  model  parameters  to  obtain  a 
realistic  picture.  (All  model  parameters  should  be  positive.)  Eight  parameters 
are  considered  as  variables  namely  b,  c,  d,  cr,  kr,  tj  and  t^  Model  parameter 
a  is  obtained  from  the  following  condition 

at 1 1 2  ♦  bt*  ♦  i  1  (10) 

Following  constraints  are  imposed  for  the  model  . 

0  <  b  <  0.3 


er  ±  c  <_  74.5 

min 

k  <  ’i  <  k 

s  .  —  s  —  s 

min  max 

k  «  k  <  k 

r  .  —  r  —  r 

min  max 

0  «.  tj  <_  1.7 
0  <_  t2  <_  1.5 

C  and  e  are  normally  set  at  0.5  and  S.0  respectively  but  Could 
min 

be  changed  in  using  input  data,  k  ,  kf  ,  k^  ,  and  k  are  given  as 

min  -min  5  max  ‘max 

data  for  the  problem.  Successful  optimization  with  least  amount  of  computer 
time  involves  selecting  a  narrow  range  for  k$  and  kr  while  keeping  it  within 
the  expected  range. 

As  the  problem  is  highly  nonlinear,  mostly  it  hits  a  local  minima,  instead 
of  true  global  minimum.  To  avoid  this  situation  the  computer  program  should 
be  run  several  times  with  different  initial  values.  Program  "Opti"  allows  to 
change  the  starting  points  and  number  of  runs  required.  It  changes  the  starting 
value  of  c,  d,  C.,  tj  and  t2  between  the  specified  values  given  as  input  dat3. 

Even  with  several  runs  of  different  starting  points  it  is  extremely  hard 
to  hit  the  local  minima.  Therefore  one  should,  select  a  local  minima  which 
gives  close  values  of  inter-  and  intra-cluster  void  ratios  for  vertical  and  horizontal 
measurements. 

This  constraint  for  selecting  local  minima  as  the  optimum  point  is  quite 
justifiable  as  irrespective  of  the  direction  of  electrical  measurement,  one  should 
get  the  same  cluster  void  ratio  even  in  highly  nonlinear  problems  like  this. 
2.3  Formation  Factor  and  Shape  Factor 

From  the  low  frequency  (1  KHZ)  conductivity  measurement*,  the  formation 
factor  F  is  defined  as  the  ratio  of  the  pore  fluid  conductivity,  ’  ,  to  the  soil 
sample  conductivity,  a,  as  follows* 
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The  formation  factor  has  been  shown  to  re.’ate  to  the  porosity  and  anisotropy 
of  soil  particles,  Arulanandan  et  ai.  (1979). 

The  formation  factor  was  shown  to  be  a  tensorial  parameter  with  r^nsorial 
components  related  to  the  microstructural  features  in  soils.  Oaf  alias  et  al.  (1979). 
The  average  formation  factor  F  and  the  anisotropy  index.  A,  are  defined,  for 
a  transversely  isotropic  soil,  as  follows: 


F  =  (Fy  ♦  2Fh)/3  (29) 

A2  =  Fy/FH  (30) 


where: 


formation  factor  in  the  vertical  direction 
formation  factor  in  the  norizontal  direction 


An  integration  technique  proposed  by  Bruggeman  (1935)  was  used  by 


Oaf  alias  and  Arulanandan  (1973)  to  derive  an  expression  for  average  formation 


factor,  F,  as  a  function  of  porosity,  n,  and  average  shape  factor,  T,  as 

F  =  n~*  (31) 

The  average  shape  factor,  I,  is  the  negative  slope  of  the  log  F  -  log  n  plot. 
It  is  this  first  invariant  of  the  second  order  shape  factor  tensor  f  and  it  relates 
the  electric  fields  inside  and  outside  the  soil  particles.  It  has  been  shown  both 
theoretically  and  experimentally  that  the  shape  factor  is  direction  dependent 
and  depends  on  porosity,  gradation  and  particl is*  shape  and  orientation, 
Kutter  (1978),  Arulanandan  et  al.  (1979),  Dafaiias  et  al.  (1979),  and  Arulmoli 
(1980).  Since  the  average  formation  factor  is  independent  of  orientation  of 


particles,  the  average  shape  factor,  for  a  given  soil,  is  expected  to  be  a  function 
of  porosity  and  the  shape  of  particles. 

The  electrical  parameters  F,  A,  f  of  soil  deposits  are,  governed  by  the 
grain  and  aggregate  characteristics  of  the  particles. 

3.  Experimental  Procedure  arid  Results 

A  detailed  experimental  study  was  conducted  to  investigate  the  variation 
of  stress  state  with  the  electrical  parameters  for  several  soils.  For  this  purpose, 
it  was  necessary  to  develope  a  series  of  identical  cells  to  measure  the  stress 
state  and  electrical  parameters  simultaneously. 

Four  identical  circular  cells  of  internal  diameter  2"  and  height  5"  (Figures  5 
and  6)  were  fabricated  out  of  teflon  tubes  of  ft”  thickness.  Two  of  the  cells 
were  used  to  measure  electrical  parameters  in  the  vertical  and  horizontal 
directions  and  the  other  two  to  measure  total  lateral  pressure  and  pare  water 
pressure. 

The  cell  used  for  the  determination  of  electrical  properties  in  the  vertical 
direction  consists  of  two  perforated  brass  electrodes  placed  at  top  and  bottom 
of  the  sample.  The  other  cell,  used  for  the  determination  of  electrical  properties 
in  the  horizontal  direction  consists  of  two  curved  brass  electrodes  of  radius  of 
curvature  1",  length  2"  and  height,  IK".  These  electrodes  were  embedded  in  the 
teflon  walls  on  diametrically  opposite  ends.  The  electrodes  in  the  cells  were 
connected  to  the  impedance  analyzer  by  a  6"  long  coaxial  transmission  cable. 

Since  the  impedance  of  the  transmission  line  and  cell  was  very  significant 
in  the  radio  frequency  range,  its  effect  had  to  be  elminiated  by  calibration. 
The  calibration  procedure  used  by  Anandarajah  (1982)  was  adopted  using  water 
of  known  conductivity  and  dielectric  constant  of  79.  Hewlett-Packard  impedance 
analyzer  (Figure  7),  which  is  designed  for  precise  impedance  measurements  in 
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the  frequency  range  of  1  MHZ  -  1000  MHZ,  was  used  for  the  measurements  in 
the  radio  frequency  range.  A  low  frequency  impedance  comparator  (Figure  8) 
was  used  for  the  determination  of  sample  conductivity  and  consequently  the 
formation  factor  -  at  a  frequency  of  100  KHZ.  The  cell  was  calibrated  to 
obtain  the  ceil  constant  using  the  low  frequency  impedance  comparator  along 
with  water  of  known  conductivity. 

Water  conductivity  and  pore  fluid  extract  conductivity  were  measured 
using  Beckmann  conductivity  meter  (Figure  9). 

Pressure  transducer  model  PG-003-5C  was  fitted  to  one  of  the  cells  for 
the  purpose  of  measuring  the  total  lateral  stress  and  it  was  in  turn  connected 
to  a  Daytronic  digital  strain  indicator  (Figure  10).  As  the  surface  of  the 
transducer  was  not  curved,  a  curved  pad  was  made  using  silicon  glue,  and  the 
transducer  was  calibrated  by  applying  hydraulic  pressure  to  the  cell.  This 
calibration  accounts  for  the  effects  of  the  pad  and  the  strains  induced  by  fixing 
the  transducer  to  the  cell. 

Pore  ■  pressure  transducer  model  GT-300  was  fitted  for  the  purpose  of 
measuring  the  pore  pressure  variation  during  consolidation.  A  small  porous  stone 
was  inserted  between  the  cell  and  the  transducer  for  the  measurement  of  pore 
water  pressure.  The  transducer  was  con  c  cted  to  the  Budd  strain  indicator 
(Figure  11).  The  transducer  was  calibrates  -  before  fixing  -  using  air  pressure 
in  the  working  pressure  range. 

The  soils  used  were: 

1.  Snow  Cal  (60%)  ♦  lllite  (40%) 

2.  Snow  Cal  (95%)  ♦  Bentonite  (5%) 

3.  Yolo  Loam 

4.  Marysville  Red  Soil 


6. 


Kadline  -  M.P. 


7.  Snow  Cal  (70%)  +  Bentonite  (30%) 

and 

Mixed  soils,  detailed  proportions  are  given  in  Table  I. 

Each  soil  except  mixed  soils  was  sieved  through  a  #200  mesh  sieve. 
Sufficient  amount  of  water  was  added  to  the  fines  to  bring  the  water  content 
to  about  3-4  times  the  liquid  limit  of  the  soil.  This  procedure  ensures  that  the 
sample  would  behave  very  similar  to  a  sample  sedimented  from  a  dilute  suspension 
Olson  (19$2).  The  slurry  was  mixed  for  6  hours  in  a  mechanical  mixer  and  kept 
inside  a  vibrator  under  a  vacuum  pressure  for  24  hours  to  ensure  complete 
saturation.  Then  the  sample  was  kept  for  2-4  weeks  for  it  to  attain  equilibrium 
with  its  pore  fluid.  The  samples  of  Yolo  Loam,  Marysville  Red  and  Kaoline  -  M.P. 
were  prepared  with  Q.05N  salt  solution. 

The  mixed  six  soils  were  prepared  by  mixing  certain  proportions  of  sand, 
silt  and  clay.  Soil  samples  were  mixed  with  a  O.05N  solution  (2.925  gm  of  NaCl 
in  10C0  cc  of  distilled  water),  and  left  for  about  twenty  hours  to  achieve  a 
state  of  equilibrium. 

The  electrical  properties  -  R  and  x  (resistance  and  reactance)  -  both  in 
horizontal  and  vertical  directions,  lateral  stress  and  exceiss  pore  water  pressure 
were  measured  at  each  stage  of  consolidation. 

The  electrical  dispersion  data  were  obtained  by  reducing  the  impedance 
analyzer  data  to  take  account  of  lino  impedance. 

The  dispersion  data  were  plotted  and  joined  by  a  smooth  curve  to  eliminate 
any  calibration  and  measuring  errors.  Subsequently  data  points  extracted  from 
those  curves  were  used  for  the  optimization  program  to  get  the  three  element 
model  parameters  and  the  intra-  and  inter -cluster  void  ratios. 


The  electrical  dispersion  curves  for  mixed  soils  at  different  water  contents 
are  presented  in  Figures  12-35  for  loading  and  unloading  condition  in  both  vertical 
and  horizontal  directions.  The  optimized  model  parameters  using  the  experimental 
dispersion  results  corresponding  to  different  water  contents  are  presented  in 
Tables  2-13.  The  formation  factor,  F,  was  plotted  against  the  total  porosity, 
n,  and  is  presented  in  Figures  36-48. 

4.  Electrical  Properties  in  Relation  to  Mechanical  Properties 

It  has  been  shown  that  the  principal  factors  influencing  the  dielectric 
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dispersion  of  fine  grained  soils  in  the  radio  frequency  range  (10  -  10  HZ)  are 
the  compositional  properties  of  the  different  phases  and  the  heterogeneous  nature 
of  tne  system.  The  engineering  properties  of  fine  grained  soils  are  also  controlled 
by  these  factors.  This  mutual  dependency  is  used  as  a  basis  to  correlate  the 
electrical  properties  to  mechanical  properties. 

The  purpose  of  this  research  is  to  establish  possible  correlations  between 
those  parameters  which  are  needed  for  the  prediction  of  stress-strain  behavior 
of  fine  grained  soils  and  the  appropriate  electrical  parameters.  It  will  then  be 
possible  to  determine  the  stress-strain  behavior  of  fine  grained  soils  by  estimating 
these  mechanical  properties  from  the  above  correlation  and  using  them  in  a 
suitable  constitutive  model.  This  approach  facilitates  the  determination  of  the 
stress-strain  behavior  of  soils  by  nondestructive  in  situ  measurements  since  the 
electrical  properties  can  be  determined  in  situ  without  disturbing  the  structure 
of  the  soil. 

Any  constitutive  model  for  the  prediction  of  stress-strain  behavior  of  fine 
grained  soils  would  require  the  knowledge  of  the  slope  of  the  consolidation  line, 
X,  the  slope  of  the  swelling  line,  «,  the  slope  of  the  critical  state  line,  M,  the 
initial  void  ratio,  eQ,  and  the  overconsolidation  ratio  (along  with  other  model 
parameters  appropriate  for  the  model  under  consideration).  These  parameters 
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are  defined  in  Figure  51.  The  water  content  of  the  fine  grained  soils  can  be 
determined  by  conventional  procedures  and  the  initial  void  ratio,  eQ,  could  be 
evaluated.  The  correlations  between  the  electrical  properties  and  the  mechanical 
properties  X,  <  and  M  for  normally  consolidated  soils  would  be  established  in 
the  following  sections. 

*.l  Slope  oi  Isotropic  Consolidation  Line,  X 

Bolt  (1956)  attempted  to  predict  the  compressibility  characteristics  of 
clays  based  on  the  concept  of  osmotic  pressure  using  Gouy-Chapman  diffuse 
double  layer  theory  and  Van't  Hoff’s  theories  of  parallel  platy  particles. 
Experimental  compression  characterics  of  Na-montmorillonite  and  Na-illite  were 
found  to  be  close  to  the  predicted  relationships.  The  theory,  however,  was 
found  to  be  valid  only  for  clays  exhibiting  very  strong  colloidal  properties  such 
as  montmorillonite.  Derivations  irom  Bolt's  findings  have  been  reported  by 
Mitchell  (1960)  and  Olson  and  Mitronovas  (1962)  and  are  ascribed  mainly  to 
particle  orientation.  Quigley  and  Thompson  (1966)  have  observed  fabric  changes 
in  natural  Leda  clay  during  consolidation  using  X-ray  diffraction  methods.  It 
has  been  shown  by  Rosenquist  (1958)  that  the  compressibility  of  clay  is  dependent 
on  the  type  as  well  as  the  valence  and  concentration  of  ions  adsorbed  on  the 
surface  of  the  clay  particles.  Further,  Olson  and  Mesri  (1970),  have  concluded 
that  both  mechanical  and  physico-chemical  factors  influence  the  compressibility 
of  soils  in  general*  although  one  or  the  other  may  dominate  depending  on  the 
soil  type. 

The  factors  influencing  the  magnitude  of  dielectric  dispersion,  A eQ,  were 
investigated  in  detail  by  Arulanandan  et  ai.  (1973).  It  has  been  found  that  AeQ 
is  significantly  influenced  by  type  and  amount  of  clay  mineral.  The  values  of 
Aeo  were  shown  to  increase  in  the  sequence  kaolinte  <  illite  <  montmorillonite. 
The  compression  index  of  these  soil  also  increases  in  this  sequence.  The 


magnitude  of  dispersion  decreases  with  an  increase  in  percentage  of  sand  in 
sand-clay  mixtures  Arulanandan  et  al.  (1973);  so  does  the  compression  index  as 
it  is  widely  known.  Olson  et  ai.  (1970)  have  shown  that  the  compression  index 
of  kaolinite  is  decreased  when  the  electrolyte  concentration  is  increased  from 
0.0001  N  sodium  to  1.0  N  sodium  and  Arulanandan  et  al.  (1973)  have  shown  that 
AeQ  also  decreased  With  increasing  electrolyte  concentration. 

The  preceding  discussion  suggests  that  the  factors  influencing  the 
mechanisms  controlling  the  compression  of  clays  and  magnitude  of  dielectric 
dispersion  are  the  same.  Based  on  this  mutual  dependency,  AeQ  has  been 
correlated  with  X.  Shading,  3.R.  (1972)  has  shown  that  there  is  a  linear 
relationship  between  AeQ  and  X  for  natural  clays  as  shown  in  Fig.  52.  Further 
results  confirm  the  general  validity  of  this  relationshi[  between  Agq  and  X. 
*.2  Slope  of  Isotropic  Swelling  Line,  ic 

The  swelling  characteristic  of  saturated  clays  due  to  the  removal  of 
external,  load  has  been  investigated  by  many,  either  by  mechanical  models  such 
as  the  one  used  by  Terzaghi  (1929)  where  swelling  is  assumed  to  result  from 
elastic  rebound  of  bent  particles  or  by  physico-chemical  models  such  as  the  one 
used  by  Bolt  (1956)  where  osmotic  repulsive  forces  are  assumed  to  be  responsible 
for  swelling.  Although  it  has  been  possible  to  explain  the  mechanism  controlling 
swelling  characteristics  by  the  above  concepts  it  was  not  very  successful  owing, 
to  the  complicated  structural  arrangements  of  particles  in  clays. 

The  concept  of  clusters  in  fine  grained  soils  (Michael  et  al.  1954  and 
Quirk,  1959)  has  beer'  utilized  by  Olsen  (1961)  in  his  study  of  hydraulic  flow 
through  saturated  clays  and  he  concluded  that  the  discrepancy  between  the 
measured  permeability  and  <  he  one  predicted  by  Kozeny-Carman  equation  in 
clays  is  mainly  due  to  unequal  pore  sizes  due  to  grouping  of  clay  particles  in 
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clusters.  The  existence  of  primary  particles  aggregation  has  been  observed  by 
many  using  electron  microscope,  Quigley  et  al.  (1966). 

Fig.  o3  shows  variation  of  intra-  (e^)  and  inter-cluster  )  void  ratios 
with  total  void  ratio  evaluated  using  electrical  dispersion  data  for  Snow  Cal 
(95%)  +  Montmorillonite  (5%).  The  results  corresponding  to  measurements  made 
in  the  vertical  and  horizontal  directions  are  identical  and  are  very  similar  to 
those  predicted  by  Olsen  (1961). 

Fig.  54  shows  for  mixed  soils  and  Fig.  55  for  clays  the  variation  of 
intra-  (ej)  and  inter-cluster  (e^)  void  ratios  with  total  void,  ratio  evaluated  using 
electrical  dispersion  data. 

It  has  been  shown  that  the  swelling  of  fine  grained  soils  is  caused  by 
swelling  .of  clusters  Smith  and  Arulanandan  (1981)  and  the  decrease  in  inter-cluster 
pores  during  compression  is  irreversible  Meegoda  (1983). 

If  the  ratio  of  intra-cluster  to  total  void  ratio  is  large  for  a  given  soil, 
the  elastic  compression  due  to  an  increase  in  the  external  load  would  be  high 
and  consequently  swelling  would  also  be  high  when  the  load  is  removed.  Assuming 
this  mechanism  of  swelling,  the  ratio,  e^/e  if  corrected  with  <  as  shown  in 
Fig.  56. 

4.3  Slope  of  Critical  State  Line  M  of  the  Bounding  Surface  Theory 

When  a  soil  element  is  sheared  under  drained  or  undrained  conditions, 
experimental  results  indicate  that  the  soil  element  fails  when  the  stress,  path 
roaches  the  critical  state  line  independent  of  the  initial  stress  state  of  the  soil 
element  (Schofield  et  al.,  1968).  At  failure,  the  void  ratio,  e,  and  the  effective 
mean  normal  pressure,  p,  lie  on  a  unique  line,  referred  to  as  a  critical  state 
line  on  the  e-p  space.  This  concept  is  widely  known  as  the  critical  state  concept 
(Roscoe  et  al.,  1968).  The  slope  of  the  critical  state  line,  on  the  p-q  space, 
therefore  represents  the  ultimate  shear  strength  of  soils. 
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Lambe  (1960)  has  discussed  the  factors  controlling  ultimate  shear  resistance 
of  fine  grained  soil  which  is  considered  to  be  due  to  friction  and  interference 
between  particles.  These  components  of  ultimate  shear  resistance  and  hence  M 
would  depend  on  many  factors  such  as  particle  size,  shape,  surface  texture  and 
the  structure  of  thj  soils  resulting  from  the  attractive  and'  repulsive  forces 
between  the  adjacent  clay  particles.  In  fine  grained  soils,  the  shape  factor 
would  reflect  in  addition  to  shape  of  the  particles  the  physico-chemical  interaction 
between  particles.  In  other  words,  the  shape  factor,  f,  is  a  function  of  the 
shape,  composition  and  arrangement  of  particles. 

Soils  which  exhibit  higher  anisotropy  were  found  to  have  lower  strength 
when  measurements  are  made  in  the  direction  of  the  major  prinicipal  axis. 
Based  on  these  considerations,  an  attempt  is  made  here  to  correlate  M  with  an 
electrical  index  defined  as  a  function  of  A  and  T. 

The  values  of  M  are  obtained  from  the  normally  consolidated  undrained 
test  results.  A  computer  program  developed  by  Herrmann  et  al.  (1980)  for  the 
calibration  of  bounding  surface  theory  (i.e.  evaluation  of  model  parameters  by 
matching  theoretical  and  experimental  stress-strain  relationships)  was  used  to 
obtain  the  values  of  M.  A  direct  assessment  from  the  failure  value  of  q/p 
could  only  be .  approximate  on  account  of  the  inaccuracies  in  measurement  of 
the  stress  parameters  at  large  strain  and  would  underestimate  M  because  failure 
intervenes  before  the  critical  state  line. 

The  correlation  between  M  and  \cfl  based  on  results  corresponding  to 
four  different  clays  and  mixed  soils  tested  is  shown  in  Fig.  57.  A  reasonable 
nonlinear  correlation  between  M  and  A^/f  is  evident. 

4.*  Correlation  of  Earth  Pressure  at  Rest,  wjth  Electrical  Parameters 
If  the  material  to  be  tested  is  purely  elastic,  Ko  will  depend  on  elastic 
properties  (Poisson's  ratio  v)  as  given  by  equation  (28),  On  the  other  hand,  if 
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Only  isotropic  materials  could  be  modelled  by  the  above  function.  Real 
soils  under  Kq  -  loading  are  anisotropic  due  to  the  presence  of  platy  type  of 
particles.  As  the  degree  of  anisotropy  increases  lateral  stresses  under  zero 
lateral  strains  decreases  since  low  lateral  stresses  are  associated  with  parallel 
arrangement  of  particles.  Therefore  mechanical  anisotropy  has  a  marked  effect 
on  Kq  in  addition  to  friction  angle  or  the  strength  of  the  material. 

Mechanical  anisotropy  could  be  quantified  by  the  electrical  anisotropy 
parameter  A.  Particle  friction  could  be  quantified  by  the  average  shape  factor 
f  as  the  shape  factor  depends  on  the  particle  roughness  Arulanandan  et  al. 
(1979). 

The  variation  of  Kq  with  consolidation  pressure  is  given  in  Figure  58. 

This  shows  that  after  a  pressure  of  around  10  psi  Kq  value  approaches  a  constant 

value  and  remains  around  that  value  on  further  consolidation. 

These  constant  values  of  K  for  each  soil  tested  are  used  to  correlate 

o 
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with  electrical  index  A  f  and  are  given  in  Figure  5?  (Meegoda  (1983)).  When 

there  is  no  friction  or  anisotropy,  is  assumed  as  1  and  the  lower  limit  of 

Kq  for  fine  grained  soils  is  assumed  as  0.38  in  deriving  the  correlation  between 

4 

K  and  electrical  index  A  i. 
o 

4 

Figure  60  shows  how  Kq  varies  with  the  electrical  index  A  f  for  different 
consolidation  pressures  for  Snow  Cal  and  5%  Montmoril Ionite.  This  result  shows 
that  the  proposed  correlation  is  justifiable  for  consolidation  pressures  greater 
than  10  psi. 

3.  In  Situ  Prediction  of  Stress-Strain  Behavior  of  Cohesive  Soils 

Application  of  the  theory  of  plasticity  has  gained  wide  popularity  over 
the  past  few  years  for  the  analytical  prediction  of  stress-strain  relationship  of 
soils.  The  most  notable  and  earliest  of  the  constitutive  models  developed  using 
the  theory  of  plasticity  is  the  Cam-clay  theory  developed  by  the  Cambridge 
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group  (Schofield  et  alM  1968).  The  Cam-ciay  model  can  represent  the  strain 
hardening  behavior  depending  on  the  state  of  the  soil  with  respect  to  the  critical 
state  line.  The  verification  of  the  Cam-clay  model  under  simple  shear  and 
triaxial  stress  states  has  been  carried  out  satisfactorily  on  soils  that  are  "wet" 
of  critical  state  (Roscoe  et  al.,  1968).  However,  this  model  has  not  been 
demonstrated  to  be  representative  of  the  strain  softening  behavior  of  the  soils. 
Yield  surface,  along  with  the  normality  rule,  gives  rise  to  prediction  of  unrealistic 
behaviors  of  isotropic  loading.  This  is  due  to  the  fact  that  the  yield  surface 
derived  as  the  result  of  assumptions  made  in  Cam-clay  theory  does  not  intersect 
the  p-axis  at  the  right  angle. 

Roscoe  et  al.  (1968)  modified  the  theory  by  introducing  an  ellipse  for  the 
yield  surface  and  using  normality  condition.  The  size  of  this  elliptic  yield 
surface  is  completely  defined  by  the  initial  isotropic  consolidation  pressure,  pQ, 
and  the  slope  of  the  critical  state  11.-10,  M.  The  value  of  R,  which  is  defined 
as  the  ratio  of  the  mean  normal  pressure,  pQ,  at  which  the  yield  surface  intersects 
the  p-axis  and  the  mean  normal  pressure,  Pj,  at  which  the  yield  surface  intersects 
the  critical  state  line  (refer  to  Figure  61),  is  assumed  to  be  2.0,  a  constant 
value.  It  is,  however,  desirable  to  be  able  to  use  different  values  of  R  in  the 
theory  in  order  to  be  able  to  predict  the  stress-strain  behavior  of  a  broad  range 
of  soil  types. 

Numerous  theories  have  since  been  developed  to  describe  the  stress-strain 
behavior  of  soils  under  a  general  loading  condition.  Bounding  surface  plasticity 
model  (Oafalias  et  al.,  1980),  cap  mode!  (Di  Vtaggio  et  al.,  1971),  endochronic 
theory  (Vaianis  et  al.,  1971),  and  model  developed  by  Prevost  (1978)  are  some 
of  them. 

The  bounding  surface  plasticity  model  has  been  developed  to  describe  a 
generalized  stress-strain  behavior  of  fine  grained  soils  under  a  variety  of 
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conditions  such  as  normally  consolidated  and  over-consolidated,  strain  hardening 
and  strain  softening,  monotonic  and  cyclic,  compression  and  extension,  loading 
paths  in  the  tension  zone  (i.e.,  the  effective  mean  normal  stress  is  tension), 
etc.  The  verification  of  this  theory  under  some  of  the  above  conditions  has 
been  carried  out  satisfactorily.  The  value  of  R  can  be  varied  as  required  to 
predict  a  given  experimental  stress-strain  relationship,  lot  like  in  the  case  of 
critical  state  theory  developed  by  Roscoe  et  al.  (1%8).  When  R  -  2,  these  two 
theories  predict  identical  stress-strain  behavior  in  the  case  of  normally 
consolidated  soils. 

).l  Bounding  Surface  Plasticity  Model 

The  theory  of  plasticity,  along  with  the  bounding  s>irface  concept  and  the 
critical  state  concept,  was  used  to  develop  a  bounding  surface  plasticity  model 
for  the  description  of  the  stress-strain  behavior  of  fine  grained  soils. 

The  bounding  surface  concept  is  that  the  plastic  modulus  of  the  material 
at  a  given  stress  state  can  be  determined  knowing  the  stress  state,  a  distance 
from  a  point  on  the  bounding  surface  (defined  by  means  of  a  suitable  mapping 
rule)  to  the  point  representing  the  current  stress  state,  a  plastic  internal  variable 
to  account  for  the  effect  of  past  loading  history  and  the  plastic  modulus 
corresponding  to  the  projected  point  on  the  bounding  surface.  The  bounding 
surface  and  a  radial  mapping  rule  are  illustrated  schematically  in  Figure  62 
where  the  bounding  surface  is  shown  in  the  space  of  Ij  (first  effective  stress 
invariant)  and  ^  (second  stress  invariant); 

The  plastic  modulus  at  point  B  can  be  found  by  the  consistency  condition 
knowing  the  functional  form  of  the  bounding  surface.  The  plastic  internal 
variable  chosen  in  developing  the  theory  for  soils  is  the  integral  of  the  plastic 
void  ratio,  e^.  The  plastic  modulus  at  the  current  state,  A,  can  be  expressed 
as  follows: 
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where  K  plastic  modulus  at  A,  =  plastic  modulus  at  B,  =  stress  tensor, 
ep  =  plastic  void  ratio,  SQ  =  reference  stress,  H  =  positive  "shape"  hardening 
function  of  the  state  of  stress. 

Knowing  the  plastic  modulus,  the  plastic  strain  increment  due  to  a  given 
stress  increment  can  be  determined  through  the  plastic  constitutive  relation 
using  a  suitable  flow  rule  as  follows: 


e. .  =  <L>n. .  (32) 

‘1  iJ 

L  =  K  °kl  nkl  (33) 

where  e..  =  plastic  strain  increment  tensor,  L  =  loading  function,  o..  =  stress 

increment  tensor  and  r»..  =  direction  of  the  plastic  strain  increment  vector. 

Assuming  normality  rule,  n.j  is  the  unit  normal  to  bounding  surface  at  the 

projection  point  B  (Figure  62).  The  behavior  of  the  material  is  elastic  in 
unloading  and  elasto-plastic  in  loading. 

One  of  the  desirable  features  of  this  approach  is  that  the  plastic 
deformation  takes  place  within  the  bounding  surface  unlike  in  the  classical 
plasticity  theory  where  the  behavior  is  fully  elastic  within  the  yield  surface. 
This  feature  of  the  model  yields  realistic  predictions  for  overconsolidated  soils 
and  under  cyclic  loading  conditions,  The  loading  function  L  (Equation  33)  includes 
K  to  account  for  loading  during  the  unstable  behavior  since  nk(  <  0  when 
K  <  0  and  consequently  L  >  0. 

The  detail  description  of  the  theory  could  be  found  in  many  references 
(e.g.,  Daf alias  et  af.  1980).  The  bounding  Surface  used  for  soils  is  shown  in  the 
triaxial  space  as  in  Figure  63. 
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The  undrained  stress  path  of  a  normally  consolidated  soil  earn  be  expressed 
by  a  closed  form  equation  as  follows: 

-3  =  [J  (JMrir  ♦  (i  - .  (f-)2]’4  04) 

R-l  K  p_  K.  P  P 

ro  ro  o  ro 

where  pQ  =  isotropic  preconsolidation  pressure,  p,  =  mean  normal  at  peak  value 
of  q  on  the  yield  line,  R  =  material  constant  =  pQ/pj. 

Since  qy  =  Mpu  at  failure,  the  undrained  strength  of  a  normally  consolidated 
soil  can  be  expressed  as  follows: 

X-ic 

%  *  M  P0  R  ’  (35) 


The  information  required  for  the  prediction  of  the  stress-strain  behavior 


of  normally  consolidated  soils  is  X,  ic,  M,  R,  eQ  and  pQ. 

5.2  Application  to  In  Situ  Prediction  of  Stress-Strain  Behavior 

The  model  parameters  defining  the  bounding  surface  plasticity  model  can 
be  obtained  using  a  set  of  triaxial  test  data  by  a  trial  and  error  procedure,  i.e., 
the  model  parameters  should  be  chosen  so  that  the  predicted  and  measured 
stress-strain  relationships  are  close  to  each  other.  A  computer  program  i  was 
developed  Herrmann  et  al.  (1982)  as  described  in  the  earlier  section  incorporating 
the  bounding  surface  theory  in  order  to  predict  the  stresses  and  strains.  This 
program  was  used  to  obtain  the  value  of  M  which  predicts  the  relationships  of 
q  vs.  c j  and  q  vs.  p  (ej  =  axial  strain)  as  close  to  the  corresponding  experimental 
relationships  as  possible.  The  values  of  "M"  obtained  in  this  manner  are  correlated 
with  the  electrical  index,  A2/f  as  shown  in  Figure  57. 

A  reasonably  good,  agreement  between  the  experimental  and  theoretical 
predictions  of  stress-strain  behavior  of  soils  is  seen  from  the  results  presented 
in  Figures  64-76.  It  is  therefore  possible  to  predict  the  stress-strain  behavior 
of  normally  consolidated  fine  grained  soils  using  bounding  surface  plasticity  model 
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if  model  parameters  can  be  determined  by  in  situ  tests.  In  situ  model  parameters 
could  be  predicted  from  in  situ  electrical  measurements.  Therefore  it  can  be 
concluded  tl  at  if  in  situ  electrical  measurements  (i.e.,  apparent  dielectric 
constants  and  conductivity  in  the  radio  frequency  range  along  with  conductivity 
of  the  solutio  extract)  are  made  in  situ,  stress-strain  behavior  and  in  situ  stress 
state  of  normally  consolidated  fine  grained  soils  could  be  predicted. 

6.  Future  Work 

It  is  proposed  to  investigate  the  possibility  of  the  prediction  of  initial  in 
situ  stress  state  using  the  electrical  method.  The  prediction  of  in  situ  stress-strain 
behavior  of  overconsolidated  soils  requires  a  knowledge  of  O.C.R.  or 
preconsolidation  pressure  (Pc)  and  the  value  of  Kq.  It  is  proposed  to  establish 
correlations  between  the  mechanical,  parameters  Kq,  Pc  and  K  and  electrical 
parameters  in  order  to  enable  the  prediction  of  KQ,  Pc  and  K  in  situ  where 

K  =  K  ,  O.C.R.®  (a  *.4505). 

o 

In  this  study,  correlations  between  electrical  and  mechanical  properties 
have  been  established  to  enable  the  prediction  of  stress-strain  relationship  of 
normally  consolidated  soils.  It  is  proposed  now  to  verify  these  predictions. 
Electrical  properties  will  be  used  to  predict  \,  k  and  M  of  normal  consolidated 
soils.  These  predicted  mechanical  properties  will  be  used  in  the  bounding  surface 
plasticity  model  to  predict  the  stress-strain  relationships.  These  relationships 
will  be  compared  with  the  measured  values.  Sufficient  number  of  soils  will  be 
tested  to  establish  the  sensitivity  of  the  predicted  parameters  X,  k  and  M  to 
the  variations  between  the  measured  and  predicted  stress-strain  relationships. 
It  is  also  proposed  to  develop  a  relationship  between  mechanical  and  electrical 
anisotropy. 
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FIG.  1  Variation  of  Dielectric  Constant  and  Conductivity  as  a  Function  of  Frequency 


Variation  of  Dielectric  Constant  and  Conductivity  as  a  Function  of  Frequency  for  Cohesive  Soils 
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Figure  ^  Three  Element  Electrical  Model 
Each  Zone  is  Represented  by  a  Parallel  Circuit  of 
Resistor  and  Capacitor .  Impedance  of  Each  Zone  is 
Determined  by  its  Dimension  (a.  b,  c.d)  and  Specific 
Conductivity  and  Dielectric  Constant  of  Material  Forming 
the  Zone  (k5.  k,.  «f) . 


Electrical  measuring  celb  and  consolidation  frames  for  the  determination  of  electrical  properties 
horizontal  and.  vertical  directions  at  different  water  contents 


Cells  and  consolidation  frames  for  the  determination  of  pore  water  pressure  and  lateral  stress  at  different 
water  content 


The  Impedance  Comparator  for  the  Measurment  of  Electrical 
Properties  at  Low  Frequency  (1  KHZ) 
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■Figure  TO  -  The  Digital  "Daytromc"  strain  indicatar  for  the  measurement  of 
lateral  stress  at  different  water  contents  during  consolidation 


Figure  11  -  The  "Budd"  strain  indicator  for  the  measurement  of  pore  water 
pressure  at  different  water  contents  during  consolidation 


Fig.  12  :  Variation  Of  Conductivity  And  Dielect 

(Soil  #1,  Vertical  Measurement,  Loadii 
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(Soil  #1,  Horizontal  Measurement,  loading) 
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Variation  Of  Conductivity  And  Dielectric  Constant  As  A  Function  Of  Frequency 
(Soil  #2,  Vertical  Measurement,  Loading) 
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Fig.  20  •  Variation  Of  Conductivity  And  Oielectric  Constant  As  A  Function  Of  Frequency 

(Soil  13,  Vertical  Measurement,  Loading) 
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Fig. 22  Variation  Of  Conductivity  And  Dielectric  Constant  As  A  Function  Of 

(Soil  13,  Horizontal  Measurement,  Loading) 
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Variation  Of  Conductivity  And  Dielectric  Constant  As  A  Function  Of  Frequency 
(Soil  #4,  Vertical  Measurement,  loading) 
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Variation  Of  Conductivity  And  Dielectric  Constant  As  A  Function  Of  Frequency 
(Soil  #4,  Vertical  Measurement,  Unloading) 
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Variation  Of  Conductivity  And  Dielectric  Constant  As  A  Function  Of  Frequency 
(Soil  #5,  Vertical  Measurement,  Loading) 
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(Soil  #6,  Horizontal  Measurement,  loading) 


Variation  Of  Conductivity  And  Dielectric  Constant  As  A  Function  Of  Frequency 
(Soil  #6,  Vertical  Measurement,  loading) 
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Figure  42  -  The  relationship  between  the  horizontal  and  vertical  formation  factors 
and  the  total  porosity,  n  for  Snow  Cal  ♦  (30%)  Montmortllonite 
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Table  10  Optimized  Electrical  Model  Parameter*  For  Soil  #3 
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Table  11  Optimized  Electrical  Model  Bwinetw  For  Soil  #3 
(Horizontal  Measurement) 
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